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We present results from simulations of two flavor QCD thermodynamics at A^j = 4 
with domain wall fermions. In contrast to other lattice fermion formulations, domain 
wall fermions preserve the full SUL(iV/) ® S\J R{Nf) symmetry of the continuum at finite 
lattice spacing (up to terms exponentially small in an extra parameter). Just above the 
C^ ! phase transition, we find that the U^(l) symmetry is broken only by a small amount. 
Y? ! We discuss an ongoing calculation to determine the order and properties of the phase 
Q ' transition using domain wall fermions, since the global symmetries of the theory are 
expected to be important here. 

r-| ■ 1. Introduction 

Xj Domain wall fermions (DWF) [^ solve the fermion doubling problem on the lattice 

D . through the introduction of a fictitious internal flavor space. After lifting the doublers 

' with a Wilson term, the remaining light fermion respects the SVi^Nf) (g) SUij(A^j) chiral 

symmetry, up to terms exponentially small in size, Lg, of the flavor space. Since the 

/\ ■ computing cost is linear in Lg it may be possible to simulate Nj = 2 QCD on the lattice 

j^ ■ with the full flavor and chiral symmetries of the continuum close to the finite temperature 

transition using today's supercomputers. 

We summarize here the efforts of the Columbia group to study two flavor QCD thermo- 
dynamics with DWF. Some of these results have also been presented elsewhere. p|J^ The 
DWF action used in this work is as in [Q with the modifications as in |^. Some details on 
the numerical methods can be found in p|. For quenched QCD thermodynamic studies 
with DWF see [^j-Q, and with overlap fermions see |T^. For a review on DWF see |TT[ 
and references therein. 

2. The U^(l) symmetry in the deconfined phase 

At low temperatures, the anomalous breaking of the 11^(1) symmetry in QCD has sev- 
eral observable consequences. For example, it explains the large mass difference between 
the 1] and rj'. At finite temperatures, the role of the chiral anomaly is not as well under- 
stood. In Nf = 2 massless QCD, if the deconfinement transition restores Ua(1) as well 
as SUi(2) ® SUij(2) then a first order phase transition is expected. On the contrary, if 
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Ua(1) remains substantially broken in the deconfinement region, universality arguments 
suggest the possibility of a second order transition. 
Earlier attempts to address this question 
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13, 14 1 using staggered fermions could not 



produce conclusive results because because of the zero-mode shifts and classical-level 
Ua(1) symmetry breaking which plague staggered fermions. At L^ = oo massless DWF 
do not break the U^(l) symmetry at the classical level, 't Hooft first recognized that the 
Ua(1) symmetry is broken by fermion zero modes which are connected to the topology 
of the gauge fields. [^ The DWF Dirac operator can have exact zero modes |l^ and, 
in an earlier study, [1^ we demonstrated that domain wall fermions are sensitive to the 
fermionic zero modes produced by smooth topological gauge configurations for fairly small 
masses and for L, > 10. 
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Figure 1. Fit of Xn - Xs = cq + 
C2mj, vol: 16^ x 4,mo = 1.9, L^ = 
16. 



Figure 2. Xtt^Xs vs. nif, vol: 16'^ x 
4,mo = 1.9, L, = 16. 



To measure anomalous symmetry breaking, we compare the susceptibilities of the pseu- 
doscalar vr mesons and scalar 6 mesons. Their difference, Xn — Xs, is plotted versus the 
dynamical bare quark mass rrif for (3 = 5.4, 5.45 in figure ^ The lattice volume is 16^ x 4, 
Ls = 16, and Pc ~ 5.325 for Nt = 4 (see section ||). The fit ansatz Xtt ~ Xs = cq + C2m^ 
is also shown. For /3 = 5.4, cq = 0.26(6) and for (3 = 5.45, cq = 0.08(3). For both fits, 
X^/d.o.f. ~ 1. We can also conclude that Ls = 16 is large enough to suppress the chiral 
symmetry breaking terms (relative to the bare quark mass), because the data fits the 
ansatz without a linear term within our current statistics. 

It is also interesting to compare the magnitude of Xn — Xs to the magnitude of the 
susceptibilities, shown in figure 0. We see that although the differences are non-zero by a 
statistically significant amount, the scale of the difference is small. It is an open question 
as to whether the small size of the U^(l) symmetry breaking is sufficient to support a 
second order phase transition. 



3. In the transition region 



Based on earlier exploratory studies using smaller 8'^ x 4 volumes, PJ^ the transition re- 
gion was localized to 5.2 < (3^ < 5.4 for the domain wall height rriQ = 1.9. To demonstrate 
critical behavior, it is important to use large spatial volumes relative to the temporal 
extent. We studied this region on 16^ x 4 volumes with rrif = 0.02 and Lg = 24. Initial 
estimates using {qq) suggested that chiral symmetry breaking effects would be no larger 
than ~ 15% throughout the transition region. The sweep through the region is shown in 
figure ^. The gauge part of the action is the standard Wilson plaquette action, indicated 
by Ci = 0. For comparison, Nf = 2 staggered data from [18| is shown. 
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Figure 3. DWF and stag (qq), vol: 
16^ X 4, ci = 0, mo = 1.9, L, = 24. 
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Figure 4. DWF {qq), vol: 16^ x 
4, ci = -0.331, mo = 1.9, L, = 24. 



In order to set the scale, we ran an 8^ x 32 simulation at /3 = 5.325 corresponding to 
the middle of the transition region. |]19[ We found that in lattice units rrip = 1.18(3) and 
mjr = 0.654(3), yielding a critical temperature T^ = 163(4) MeV and m^ = 427(11) MeV. 
The critical temperature is in agreement with results obtained from other fermion reg- 
ulators. ||2^ The pion mass is clearly too heavy to be able to extract useful information 
regarding the order of the transition. ( By comparison, for staggered fermions near T^ for 
Nt = 4, one pion has a mass of m^ ~ 228 MeV and the other two are tti^j ~ 604 MeV, 



18 1 which might help explain the different crossover rates between the two regulators.) 



Furthermore, a more sophisticated study of residual chiral symmetry breaking effects in- 
dicate that effects of finite Lg are of the same magnitude as the bare quark mass rrif and 
suggest that Lg ~ 100 may be necessary to lower the three pion masses to the physical 
regime. ^1\ 

In an attempt to reduce the chiral symmetry breaking effects due to finite Lg, we 
improved the gauge action by adding a 1 x 2 rectangular plaquette to the standard Wilson 
gauge action, with the choice of coefficient, ci = —0.331, suggested by Iwasaki. [^ An 
initial study of the quenched hadron spectrum at couplings corresponding to the quenched 
Nt = 4 QCD transition region indicated a significant reduction in finite Lg effects. |1T9[ 
These results prompted a study of the Nf = 2 QCD transition region using Iwasaki 
improved gauge action and DWF. The results are presented in figure ^. 

Again, we set the scale simulating meson masses on an 8^ x 32 lattice at /3 = 1.9, 



corresponding to the middle of the transition region. [T^ In lattice units, we found 
rup = 1.163(21) and m^ = 0.604(3), which gives T^ = 166(3) MeV and m^ = 400(7) MeV. 
The critical temperature is in agreement with the standard Wilson gauge action results. 
However, the pion mass is not significantly lighter. 

4. Conclusions 

We measured the difference Xn — Xp just above the Nt = 4 finite temperature transition 
and found to be non-zero by a statistically significant amount. However, the magnitude 
of the difference is much smaller than the susceptibilities themselves. Thus, while the 
Ua(1) is broken in the symmetric phase it is not clear what effect this may have on the 
order of the transition. 

We studied the transition region using two different gauge actions. In both cases, the 
transition appears to be a smooth crossover, probably due to larger than expected chiral 
symmetry breaking effects of finite Lg. We then ran zero temperature simulations of meson 
masses at couplings corresponding to the middle of the crossover region to set the scale. 
In each case, the critical temperatures and pion masses we extracted are comparable to 
other lattice regulators. Alternate techniques for reducing chiral symmetry breaking at 
finite Lg are required to reduce the pion masses to physical values as simply repeating the 
study at larger Lg would take several years using current supercomputers. 

All calculations were done on the 400 Gflops QCDSP machine at Columbia University. 
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